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Abstract 

Radio-frequency El transitions between nearly degenerate, opposite parity levels of atomic dyspro- 
sium were monitored over an eight month period to search for a variation in the fine-structure constant. 
During this time period, data were taken at different points in the gravitational potential of the Sun. 
The data are fitted to the variation in the gravitational potential yielding a value of ( — 8.7 ± 6.6) x 1CF 6 
for the fit parameter k a . This value gives the current best laboratory limit. In addition, our value of 
k a combined with other experimental constraints is used to extract the first limits on k e and k q . These 
coefficients characterize the variation of m e /m p and m q /m p in a changing gravitational potential, where 
m e , m p , and m q are electron, proton, and quark masses. The results are k e = (4.9 ± 3.9) x 10~ 5 and 
k q = (6.6 ±5.2) x 10" 5 . 

According to general relativity, values of fundamental constants are independent of space and time as re- 
quired by the Einstein Equivalence Principle (EEP). Modern theories that aim to unify gravitation with the 
other forces, however, do not set this restriction on the fundamental constants; so, they may vary either 
spatially or temporally, which is in conflict with EEP pQ. 

Various recent studies have reported the results of searches for a temporal variation of the fine-structure 
constant, a = e 2 /he. These include the analyses of absorption spectra from quasars [2-5], the analyses of 
the nuclear products of the natural fission reactor at Oklo, which operated 2 x 10 9 years ago [6-9], and var- 
ious laboratory investigations involving atomic clocks, see for example Rcfs. [10-13]. In contrast to studies 
involving analyses of the processes that have occurred billions of years ago, laboratory searches are sensitive 
to present-day variation of fundamental constants. Recently, our group's experiment utilizing the El radio 
frequency (rf) transitions between nearly degenerate opposite parity levels in atomic dysprosium has yielded 
a result of a/a = (— 2.7±2.6) x 10 -15 /yr [TJj over an eight month observation period. A detailed description 
of the experimental setup and analysis was given in Refs. [14-16]. Here, we investigate a possible correlation 
between a change in a and a change in the gravitational potential. 

The energy of an atomic level can be written as 



where Eq and ao are the present-day values of the energy and a, respectively, and q contains the sensitivity 
of the level to a |17| . From Eq. (1), a change in a should result in a change in v, for a « ao, as given by 
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Similarly, the change of the rf transition frequency between levels A and B, as shown in Fig. [TJ is given by 

§{Au) = 2 (^_^1 *± , (3 ) 
h a 

where \2(qs — qA)/h\ « 1.8 x 10 15 Hz [17]. Here, we take this number to be exact, but the uncertainty is 
expected to be approximately twenty percent. 
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Figure 1: The relevant energy-level, population, and detection scheme for the dysprosium experiment. 



The transitions utilized in the experiment are the 235-MHz (J = 10 — > J = 10) transition of the isotope 162 Dy 
and the 3.1-MHz (F = 10.5 -» F = 10.5) transition of 163 Dy. The long-lived level B (lifetime r > 200 fis [18]) 
is populated via three transitions. Laser light at 833 nm and 669 nm is used for the first two transitions, 
and then the atoms spontaneously decay to level B as shown in Fig. [TJ The transition between levels A and 
B is induced by an applied rf electric field. The atoms decay from level A in two steps, and the fluorescence 
at 564 nm is monitored to detect the rf transition. 



Data were taken over an eight month period and hence at different points in the gravitational potential 
of the Sun. The gravitational potential at the Earth due to the Sun [19] is given by 



U 



-GM S 



(4) 



where G is the gravitational constant and M s is the mass of the Sun. For an elliptical orbit, the distance r 
between the Sun and Earth is 



1 



1 + 6 COS ( 



(5) 



where a is the semi-major axis of the Earth's orbit, e = yl — b 2 /a 2 « 0.0167 is the eccentricity, b is the 
semi- minor axis, and <j> is the true anomaly (see Fig. [2]). Substituting Eq. (5) in to Eq. (4), the gravitational 
potential becomes 



U = 



-GAL GAL 



:cos0 + O(e 2 ) 



(G) 



a a 

The first term in Eq. (6) is the constant part of the gravitational potential at Earth, while the second term 
is a change in gravitational potential, AU, which arises due to the eccentricity of the Earth's orbit. The 
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Figure 2: Left: The true anomaly, 0, is the angle subtended from perihelion, the point of closest approach. 
The eccentric anomaly, e. is the angle between perihelion and the position of the Earth in its orbit projected 
onto the auxiliary circle of the ellipse (the eccentricity of the ellipse is exaggerated for clarity). Right: The 
change in gravitational potential of the Sun at the Earth due to the cllipticity of the orbit. The solid line 
represents the time period during which data were taken. 



fractional change in a as a function of gravitational potential can be parametrized as |20j 

* = ijm, (7) 

a c z 
GM 

AU(t) = ecos cj)(t) . (8) 

a 

The data were taken between October 2005 to June 2006; the true anomaly is zero at the 2005 perihelion, 
January 2 [21] . The true anomaly was calculated for each data point using a two step process. The elapsed 
time in days from perihelion is related to the eccentric anomaly (see Fig. [2]) as [22] 



t = V^w :< - £_esine ' > ' ^ 

where M = M s + Me is the sum of the masses of the Sun and Earth. The eccentric anomaly is obtained by 
solving Eq. (9) and then used to calculate the cosine of the true anomaly as a function of time since 

cos e — e ,„ . 

coscj) = . (10) 

1 — e cos e 

The calculated values for cos cf) are substituted into Eq. (8) to find the gravitational potential for each 
data point. The measured frequencies for each isotope are fitted to the gravitational potential using a 
two-parameter least squares fit given by 

S(AV) -V* ^Xo^^+X! , (11) 

where v* is an arbitrary reference frequency, and Xq and X\ are the fit parameters. The parameter x\ accounts 
for the offset due to the reference frequency while the parameter x a determines the correlation between the 
change in gravitational potential and the transition frequency. If such a correlation is due to a change in a, 
the sum and difference frequencies may also be fit to the varying potential to extract additional information 
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regarding the variation of a. Since the energy difference between levels A and B is of opposite sign for the 
two isotopes and the q values are independent of the nucleus, the sum of the transition frequencies should 
be insensitive to a variation, while the difference of the transition frequencies should be twice as sensitive 
[14,15]. The data plots and the fits of the measured frequencies for each isotope as well as the sum and 
difference frequencies are shown in Fig. [3J 



235 MHz Transition in Varying Grav. Pot. 3.1 MHz Transition in Varying Grav. Pot. 
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Figure 3: Top left: The data fitted to the gravitational potential for the 235-MHz transition. Top right: The 
data fitted to the gravitational potential for the 3.1-MHz transition. The bottom plots show the difference 
and sum frequencies fitted to the change in gravitational potential. 



The parameter Xq obtained by the least-squares fit is (—0.2 ± 1.7) x 10 10 Hz for the 235-MHz transi- 
tion, (2.9 ± 1.7) x 10 10 Hz for the 3.1-MHz transition, (2.7 ± 2.4) x 10 10 Hz for the sum frequency, and 
(—3.1 ± 2.4) x 10 10 Hz for the difference frequency. 

These values can be used to calculate the constraint on the parameter k a from Eq. (7). Substituting 
Eq. (3) into Eq. (7), we get a relation similar to Eq. (11), 

6(Av) = (2 q J^y a ^± . (12) 

The fit parameter yields k a = (-1.1 ± 9.2) x 10" 6 for the 235-MHz transition, k a = (-16.3 ± 9.4) x 1CT 6 
for the 3.1-MHz transition, k a — (—8.7 ± 6.6) x 10~ 6 for the difference frequency, and (—7.5 ± 6.6) x 1CP 6 
for the sum frequency. Although the value for the sum frequency is expected to be consistent with zero, the 
one-sigma mismatch of our value is consistent with the estimated uncertainties, dominated by systematic 
effects QI]. 
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The results relating a changing a to a changing gravitational potential are summarized in Table Q] The 
results in the first row of Table [1] depend on the variation of a as well as the variation of m e /m p . The 
dependence on m e /m p is contained in k e |20j . The results in the second and third rows in Table [T] depend 
on the variation of a and the variation of m q /m P7 which is contained in k q [20| . In addition to the results 
shown in Tabic [1] our results may limit the parameter range in some specific theories, such as those in Rcfs. 
[23-25]. 

Our calculated value of k a can be combined with the results for k a + 0.17fc e and k a + 0.13fc 9 to extract 
a value of (4.9 ± 3.9) x 1CT 5 for fc e and (6.6 ± 5.2) x 1CT 5 for k q . 



Parameter 


Constraint 


Experimental Rcf. 


k a + 0.17fc e 


(-3.5 ±6) x l(r Y 


26] 


\k a + 0.13k q \ 


< 2.5 x l(T b 


m\ 


k a + 0.l3k q 


(-1 ± 17) x 1(T Y 


m 


k a 


(-8.7 ±6.6) x l(r b 


this work 


ke 


(4.9 ±3.9) x lCT b 


this work 


k q 


(6.6 ±5.2) x 10~ b 


this work 



Table 1 : A summary of results for changing fundamental constants in a varying gravitational potential based 
on the theoretical work from Ref. [2"0] . 

In this paper, we have reported the first laboratory result for k a and subsequently extracted a limit on k e and 
k q . In principle the results from Rcf. [26j can also be combined with other optical-Cs clock comparisons in 
such a way as to extract k a independent of k e . Currently, our results are limited by systematic uncertainties 
|14j . A new apparatus is under construction to address the systematic effects and, thus, possibly increase 
the sensitivity of the experiment by up to three orders of magnitude |15j . 

The authors are grateful to M. Kozlov for stimulating discussions. This work has been supported in part 
by Los Alamos National Laboratory LDRD, NSF REU supplement, and by grant RFP1-06-15 from the 
Foundational Questions Institute (fqxi.org). 
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